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Abstract 



A relation between the eigenvalues of an effective Hamilton operator and 
the poles of the S matrix is derived which holds for isolated as well as for 
overlapping resonance states. The system may be a many-particle quantum 
system with two-body forces between the constituents or it may be a quantum 
billiard without any two-body forces. Avoided crossings of discrete states as 
well as of resonance states are traced back to the existence of branch points in 
the complex plane. Under certain conditions, these branch points appear as 
double poles of the S matrix. They influence the dynamics of open as well as 
of closed quantum systems. The dynamics of the two-level system is studied 
in detail analytically as well as numerically. 



I. INTRODUCTION 

Recently, the generic properties of many-body quantum systems are studied with a re- 
newed interest. Mostly, the level distributions are compared with those following from 
random matrix ensembles. The generic properties rule, well expressed in the center 

of the spectra where the level density is high. As an example, the statistical properties of the 
shell-model states of nuclei around ^^Mg are studied a few years ago by using two-body 
forces which are obtained by fitting the low-lying states of different nuclei of the 2s — Id 
shell. In the center of the spectra, the generic properties are well expressed in spite of the 
two-body character of the forces used in the calculations. 

Another result is obtained recently in performing shell-model calculations for the same 
systems with random two-body forces. In spite of the random character of the forces, the 
regular properties of the low-lying states are described quite well in these calculations. 
Further studies proved the relevance of the results obtained and could explain in detail 
even the regular properties at the border of the spectra obtained from random two-body 
forces [|]. The spectral properties of the two-body random ensemble studied 30 years ago 
are reanalyzed 0. 

The spectra of microwave cavities are not determined by two-body forces. Nevertheless, 
the calculated spectra are similar to those from nuclear reactions ||^. They show devi- 
ations from the spectra obtained from random matrix theory as well as similarities with 
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them. Avoided level crossings play an important role. The theoretical results obtained are 
confirmed by experimental studies 

The effect of avoided level crossing (Landau-Zener effect) is known and studied theo- 
retically as well as experimentally for many years. It is a quite general property of the 
discrete states of a quantum system whose energies will never cross when there is a certain 
non-vanishing interaction between them. Instead, they avoid crossing in energy and their 
wave functions are exchanged when traced as a function of a certain tuning parameter. The 



avoided level crossings are related to the existence of exceptional points |1I0|. The relation 
between exceptional points and the well-known diabolical points is not investigated up to 
now. The relation of the latter ones to geometrical phases is studied experimentally |]1T 



as 



well as theoretically [|T2| , |l3i in a microwave resonator by deforming it cyclically. The results 



show non-trivial phase changes when degeneracies appear near to one another |T^. The in 



fluence of level crossings in the complex plane onto the spectra of atoms is studied in |TJ] . In 
this case, the crossings are called hidden crossings [T^. The relation between avoided level 
crossings and double poles of the S matrix is traced in laser-induced continuum structures 
in atoms 

Usually, it is assumed that avoided level crossings do not introduce any correlations 
between the wave functions of the states as long as the system parameter is different from 
the critical one at which the two states avoided crossing. Counter-examples have been found, 
however, in recent numerical studies of the spectra of microwave cavities P|,|l8| . These results 



coincide with the idea that avoided crossings are a mechanism of generating the random 



matrix like properties in spectra of quantum systems. 

It is the aim of the present paper to study in detail the dynamics of quantum systems 
which is caused by avoided level crossings. They are traced back to the existence of branch 
points in the complex plane where the interaction between the two states is maximum. The 
wave functions are bi-orthogonal in the whole function space without any exception. 

The paper is organized as follows. In section 2, the relation between the eigenvalues of 
an effective Hamilton operator and the poles of the S matrix is derived. This relation holds 
also in the region of overlapping resonances. It is used in section 3, where the mathematical 
properties of branch points in the complex plane and their relation to avoided level crossings 
and double poles of the S matrix are sketched by means of a two-level model. In section 
4, numerical results for states at a double pole of the S matrix as well as for discrete and 
resonance states with avoided crossing are given. The influence of branch points in the 
complex plane onto the dynamics of quantum systems is traced and shown to be large. The 
results are discussed and summarized in the last section. 



II. HAMILTON OPERATOR AND S MATRIX 

A. The wave function in the space with discrete and scattering states 

The solutions of the Schrodinger equation 

iH-E)<iI% = (1) 

in the whole function space contain contributions from discrete as well as from scattering 
states. The discrete states are embedded in the continuum of scattering states and can 
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decay. In the following, we will represent the solutions by means of the wave functions 
of the discrete and scattering states. 

Following pO| , pT| , we define two sets of wave functions by solving first the Schrodinger 
equation 

(H^' - Et) K = (2) 
for the discrete states of the closed system and secondly the Schrodinger equation 

Y^{H-' - E) = (3) 

d 

for the scattering states of the environment. Here, is the Hamilton operator for the 
closed system with discrete states and H'^'^ is that for the scattering on the potential by 
which the discrete states are defined, E is the energy of the system and the channels are 
denoted by c. 

By means of the two function sets obtained, two projection operators can be defined, 

N ^ roc 

Q = E P = T. dE (4) 

R=l c=l •^'^'^ 

with Q ■ Ce^'^ = ; P • = 0. We identify H^^ with QHQ = Hqq and H"""' with 
PHP = Hpp. From (|I]), it follows 

{Hpp - E) ■ P^t^ = -HpQ ■ gvl>g+) ; {^Hqq - E) ■ Q^f~^ = -Hqp ■ Pvl>g+) (5) 

and 

where HpQ = PHQ and Hqp = QHP. Further, 

G^^^ = P{E^+^ - HppY^P (7) 
is the Green function in the P subspace and 

H""^ = Hqq + HqpG^p^ HpQ (8) 

is an effective Hamiltonian in the function space of discrete states. 
Assuming Q + P = 1, it follows from Eq. 

^% = {P + Q) ^% = + (1 + GpHpq) ■ Qm% (9) 

where the ^'^^ on the wave functions are omitted for convenience. Using the ansatz 

Qm% = Y.BRK (10) 

R 

and Eq. (H), one gets 
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R' 

and 



n = CE + + ^R){<\^:^\^t')«'\HQP\CE) (12) 

R.R' 

where the 

un = GpHpQ^i (13) 
follow from the solutions of the Schrodinger equation 

- E) {Ce*M = {CElHpQl^t) (14) 

c' 

with source term which connects the two sets {^p} and {^e} of wave functions. With these 
coupling matrix elements 

W^n = {eE\HpQ\<^t) = {<l>t\HQp\CE) (15) 
between discrete states and scattering wave functions, it follows from (^) 



{<\H''^\<') = (<^>5il^K) + / dE' -^vr Y^W^R^h ■ (16) 

fee 

The principal value integral does not vanish, in general. 

With the eigenfunctions $r and eigenvalues £r = Ep — ^Tp of H'^^, the solution \E'|; of 
the Schrodinger equation in the whole function space of discrete and scattering states reads 

n = eE + T.^R ^ • (17) 

In order to identify the eigenvalues and eigenfunctions of H^^ with values of physical rele- 
vance, the two subspaces have to be defined in an adequate manner. When the P subspace 
contains all scattering states defined by their asymptotic behaviour and the Q subspace is 
constructed from all the wave functions of the closed system in a certain energy region (for 
details see p0|j2l[| ), the values 

W^R = rE\HpQ\^p) = mHQp\i%) (18) 

are the coupling coefficients between the resonance states and scattering wave functions, 
while the eigenvalues determine the energies Ep and widths T p of the resonance states. The 

np = ^R + Cjp = {1 + GpHpq)^p (19) 

are the wave functions of the resonance states (with ujp defined by Eq. ([T3|) when is 
replaced by '^p). 
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The Hamiltonian H*^^ is non-Hermitian since it is defined in a subspace of tlie wliole 
function space. Tlie left and right eigenf unctions, and $5?; of ^ non-Hermitian matrix 
are different from one another. For a symmetrical matrix, it follows 

(I-^jI = ($^1 and H''''\^r)=Sr\^r), (20) 



see e.g. [^|r7|,^^. Therefore, = $5* * = The eigenfunctions of if^** can be 

orthonormalized according to 

= {^},\^n') = Srr, (21) 
where $5^/ = ^r/ . Eq. ( PT| ) provides the bi-orthogonality relations 

{^r\^r) = ^{{^r\^r)) = {^R'l^R') ; Ar = {^r\^r) > 1 
= ' "^{{^rI^R'^r)) = -{^r'^rI^r) ; bS^"" = \{^r\^r'^r)\ > . (22) 

Using the orthonormality condition (0), it follows 

f^^ = -2$5{(l>^|/f<=ff|<l«)} = 27r j2iW'Ry (23) 

c=l 

for the relation between the Tr and PV]^. This relation holds also for overlapping resonance 
states. Here (H) holds, but Tr = {2tt/Ar) E \W^\'^ < 2ttJ2\W^\'^ according to (||). 



It should be underlined here that the expression ( pTD is obtained by rewriting the 
Schrodinger equation (|l[) with the only approximation P -\- Q = 1. The and as well 
as the Qr, Er, Tr depend on the energy E of the system. The energies Er = Er{E = Er), 
widths Tr = Tr{E = Er) and wave functions Qr = (Ir^E = Er) of the resonance states 
can be found by solving the corresponding fixed-point equations. Also the coupling ma- 
trix elements are complex and energy dependent functions, generally. For numerical 



examples see 



The expression (|T^ is solution of Eq. (|T|) independently of whether or not the Hamilton 
operator H contains two-body residual forces V. It is H = Hq + V e.g. in nuclear physics 
but H = Ho for quantum billiards. When H = Ho + V,it follows = {^r\V\^%) 
In the case of quantum billiards, the can be calculated by using Neumann boundary 
conditions at the place of attachment of the lead to the cavity whereas Dirichlet boundary 



conditions are used at the boundary of the cavity (see [0). In this case, the are real, 
i.e. the principal value integral in ([TBI) vanishes [|26|. The may be complex, nevertheless 



25 . 



B. The S matrix 

The S matrix is defined by the relation between the incoming and outgoing waves in the 
asymptotic region. Its general form is 

S,,, = exp{2i6,)6,,, - 2i'K{xi*\V\^!%) (24) 

where the x'e uncoupled scattering wave functions obtained from 
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$:(i/o^^' -E)xi = ^. (25) 



The are given by Eq. ([T7|). When the two subspaces are defined consistently, Eq. (pif ) 
can be written as 

Scd = - Sf] (26) 

where 

S^, = exp{2t6,)6,,, - 2z7r{xi*\VrE) (27) 
is the smooth direct reaction part and 

= 2^vr E {xi*\V\Qn) ■ ^ - (28) 

R=l ^ — + 2^ R 

is the resonance reaction part. In S^l] , V may be zero (no couphng between the channels). 

Since Eq. (|19D between the resonance states Qr and the eigenf unctions of H'^^ is 
completely analogous to the Lippman-Schwinger equation 

C^ = {1 + GpV)xe (29) 

(which describes the relation between the two scattering wave functions and Xe with and 
without channel-channel coupling, respectively), one arrives at ET 



{x'S\V\nn) = {CE\V\^R) = W'R- (30) 

When V = 0, channel-channel coupling may appear due to the coupling of the scattering 
wave functions via the Q subspace (an effective Hamilton operator in the P subspace can 
be derived analogously to the effective Hamilton operator H'^^ in the Q subspace, Eq. (^, 
see 



Using this relation, the resonance part (|28D of the S matrix reads 



SH) = 2„ ± . (31) 

R=l ^ — ^R + 2^ R 

Thus, the poles of the S matrix are the eigenvalues Sr = Er — | F/j of the Hamiltonian H*^^ 
at the energy E = Er (solutions of the fixed-point equations). The numerator of Eq. (|3lD 
contains the squares ($r)^ of the eigenfunctions of H'^^ (and not the |$i?,P) which are always 
finite in accordance with the normalization condition (^Tj). Numerical results for strongly 
overlapping resonances can be found in [E2 . 



For overlapping resonance states (i.e. W^{E)W^,{E) 7^ at the energy E for R ^ i?'), 
the W may be very different from the W. Even when the W are real, the W may be 
complex since the eigenfunctions of if ^ are complex. The coupling strength of the system 
to the continuum of channel wave functions is given by the sum of the imaginary parts of 
the diagonal matrix elements or of the eigenvalues of H'^^, 
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R R Rc Rc 



(32) 



where further Eq. (p3D is used. All redistributions taking place in the system under the 
influence of a certain parameter must obey the sum rule (p2|) . According to this rule, 
resonance trapping may appear when the resonance states overlap, 

iV _ Af _ ^ _ 

Y,VR^Y.^n : E r^^O. (33) 

_R=1 _R=1 R=M+1 

It means that, under certain conditions, N — M resonance states may decouple from the 
continuum of scattering states, i.e. they may be trapped by M states. For numerical results 
on nuclei see and on open quantum billiards see |^ . Studying the system by means of 
its coupling to the environment (described by the coupling matrix elements Wp^) will give, in 
such a case, information either on the N — M long-lived states on the background of the M 
short-lived states or on the M short-lived states with fluctuations arising from the N — M 
long-lived states. 

This behaviour induced by the imaginary part of the non-diagonal matrix elements of (|16|) 
differs from that induced by the real part. While the imaginary part causes the formation 
of structures with different time scales (as discussed above), the real part causes equilibrium 
in time (approaching of the decay widths). The first case is accompanied by an approaching 
of the states in energy {clustering of levels) while the second case is accompanied by level 
repulsion in energy. Numerical examples can be found in [l^ll^] for atoms and in ||^ for 
quantum billiards. The interplay between the real and imaginary parts of the non-diagonal 
matrix elements of (|16D characterizes the dynamics of open quantum systems. 

Eq. (|3ll) coincides formally with the standard form of the resonance part of the 5* 
matrix. It should be underlined, however, that the W^, Er, and F/j are not parameters 
(as in the standard form, see [^|), but energy dependent functions which can be calculated. 



Eq. ( |3TD holds also for overlapping resonances. Due to the energy dependence of the Tr, 
the line shape of the resonances differs from a Breit-Wigner shape, as a rule, even without 
any interferences with the direct reaction part. Some numerical examples are discussed in 



III. BRANCH POINTS IN THE COMPLEX PLANE 

Eq. (^) gives the S matrix for isolated as well as for overlapping resonance states. The 
extreme case of overlapping corresponds to the double pole of the S matrix at which the 
eigenvalues £1^2 of two resonance states are equal. Let us illustrate this case by means of 
the complex two-by-two Hamiltonian matrix 





W = ' ' - M ' ' ■ (34) 



The unperturbed energies and widths 7^ {k = 1, 2) of the two states depend on the 
parameter a to be tuned in such a manner that the two states may cross in energy (and/or 
width) when u = 0. The two states interact only via the non-diagonal matrix elements u 
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which may be complex, in general, see Eq. ([161) . In the following, we consider real uo and 7a; 
independent of a. 

The eigenvalues of Ti are 

^ ei + £2 1 



^ E,, - -r,, = ± - V(6i - e,f + 4o;2 (35) 



with i,j = 1,2 and = — | 7^ (A; = 1,2). According to Eq. (|35|) , two interacting 
discrete states (with 7^ = 0) avoid always crossing since uj and (ei — €2) are real in this 
case. Eq. (pSf ) shows also that resonance states with non-vanishing widths 7fc avoid mostly 
crossing since 

F{a,uj) = {ti-t2f + Auj^ (36) 

is different from zero for all a, as a rule. Only when F{a^uj) = at a = a'^'' (and uj = uj^"^), 
the states cross, i.e. £1 = £2- In such a case, the S matrix has a double pole, see e.g. p9| . 
It can further be seen from Eq. (P5|) that the crossing points are branch points in the 



complex plane. The branch point is determined by the values {uY and (ei — €2)^ but not by 
the signs of these values. According to Eq. (P3|), it lies in the complex plane at the point 
X = (l/2){ei(a'^'') + 62(0^"^)}. According to Eq. (pi]), it is a double pole of the S matrix. 

The eigenfunctions can be represented in the set of basic wave functions $^ of the 
unperturbed matrix corresponding to u; = 0, 

<^'.=E%*;- (37) 

In the critical region of avoided crossing, the eigenfunctions are mixed: \hii\ = \bjj\ 7^ 1 and 
bij = —bji 7^ for i 7^ j. The bij are normalized according to Eq. (|2l|) . 
The Schrodinger equation with the Hamiltonian Ti can be rewritten as 




= J2 {^k\W\^,) J2 {^k\^m)\^m) . (38) 
k=l,2 m=l,2 

Eq. ( p8| ) is a Schrodinger equation with the Hamiltonian Hq of the unperturbed system 
(corresponding to u = 0) and a source term which is related directly to the bi-orthogonality 
of the eigenfunctions of the Hamiltonian Ti. It vanishes with = ($j|<l>j) 1 and 

^ 0. Using Eq. (|2|), the rhs of Eq. (HD reads 

W\^i) = W^'[A\^i) + z5|$2)) + W^'[A\(!>2) - (39) 

with W^^ = k = 1,2. At the branch point, the complex energies are equal, 

£1^ = £2^, and 



= W,'i{A\<!>f) - tB\<^\^)) + W^i{A\^'r) + ^51$^)) (40) 
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where the relations W^^ = W"^"^ ; W^'^ = W'^^ are used. This gives 

{A + zB)\<l>'r) = {A - ^B)\<I>',^) (41) 

and finally 

^ ± i \t%d (42) 



by using the bi-orthogonality relations (^4 ) 

|($?n$^OI-oo |(«f^1$S,)|-oo (43) 

at the branch point. The relation ( ^2]) between the two wave functions at the branch point 
in the complex plane has been proven in numerical calculations for the hydrogen atom with 



a realistic Hamiltonian [r7|. Note, that the condition (21) 



($fl$?) = 5.,. (44) 

is fulfilled also at the branch point in the complex plane. This is achieved since the difference 
between two infinitely large values may be (for i 7^ j) or 1 (for i = j). Numerical examples 
for the values with z = j as well as with i j can be found in p^j23|j30| -p2| . 



Thus, the condition (^) is fulfilled in the whole function space without any exception. 
According to Eq. (|1^), the two wave functions can be exchanged by means of interferences 
between A and B, i.e. by means of the source term in the Schrodinger Eq. (|38|). When the 
wave functions are exchanged, 3fJ(6ii) = 3^(622) as well as 3fJ(6i2) = — 3f^(&2i) do not change 
their signs while 3(6ii) = $5(622) and $5(612) = —$5(621) jump between —00 and +00 at the 
branch point. 

Note that Eq. (E2|) differs from the relation = ±$2^ used in the literature for the 



definition of the exceptional point H^. With = ±$2^, neither ($^^|$^^) = 6ij nor 



= 6ij can be fulfilled. With = ±1^2^ however, the orthogonality relations 
(21) are fulfilled in the whole function space, without any exception. 

The effects arising from the source term in Eq. (^) play the decisive role in the dynamics 
of many-level quantum systems caused by avoided level crossings. This will be illustrated 
in the following section by means of numerical results. The effects appear everywhere in 
the complex plane when only 7^ 1, 7^ 0. They appear also in the function 

space of discrete states where the are real, due to the analyticity of the wave functions 
and their continuation into the function space of discrete states. 



IV. NUMERICAL RESULTS 

The numerical results obtained by diagonalizing the matrix (Q) are shown in Figs. |I| to 
^. The Ei and Fj are in units of a chosen arbitrarily, and the bij are dimensionless. In all 
cases ei = 1 — a/2, 62 = a and u = 0.05. The 7^ do not depend on the tuning parameter 
a. The relation between them is 72 = I.I71. At a = a" = 2/3, the two levels cross when 
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unperturbed (i.e. c<j = 0) and avoid crossing, as a rule, when the interaction uj is different 
from zero. Here ci = 62 = ef^^ 2 = 2/3 and 

F"(a, u) = 3ft(F"(a, u)) = Au' - (1/4) (7. - 7.)' (45) 

according to Eqs. (|35D and (|36D . The F^^{a, u) may be positive or negative. Thus, either the 
widths Ti j or the energies Eij cross freely at a^"^, but not both. The only exception occurs 
when the S matrix has a double pole at a", i.e. when 71/2 = = 1.0, 72/2 = 72^^/2 = 

1.1. Here, F'^^{a,uj) = and the two resonance states cross in spite of 7^ 0. According to 
Eq. (|35D , the double pole of the S matrix is a branch point in the complex plane. 

In Fig. |T], the energies Eij, widths Tij and wave functions bij of the two states are 
shown as a function of the parameter a in the very neighbourhood of the branch point. 
Approaching the branch point at a", |3?(6jj)| — >• 00 and |3(6jj)| — » 00. While 3?(&ij) does 
not change its sign by crossing the critical value a™, the phase of Q{bij) jumps from ± to 
=1=. The orthogonality relations (^1]) are fulfilled for all a including the critical value a™. 

Fig. 1^ shows the energies Eij and widths Tij of the two states for values of •jij just above 
and below the critical values 7^^ as well as for 'jij = 0, i.e. for discrete states. According 
to Eq. (^31), either the energy trajectories or the trajectories of the widths avoid crossing at 
the critical value a'^^ (since the condition for the appearance of a double pole of the S matrix 
is not fulfilled). It is exactly this behaviour of the trajectories which can be seen in Fig. ^ 

- When 7j > 'jf^, the widths of the two states approach each other near a'^'' but the width 
of one of the states remains always larger than the width of the other one. The two 
states cross freely in energy, and the wave functions are not exchanged after crossing 
the critical value a'^''. 

- The situation is completely different when 7j < 7". In this case, the states avoid 
crossing in energy while their widths cross freely. After crossing the critical value a'^'", 
the wave functions of the two states are exchanged. An exchange of the wave functions 
takes place also in the case of discrete states (7^ = 0). This latter result is well known 
as Landau-Zener effect. It is directly related to the branch point in the complex plane 
at a'^'' as can be seen from Fig. 0. 

The wave functions bij are shown in Fig. ^. The states are mixed (i.e. \bii\ 7^ 1 and 
bij^i 7^ 0) in all cases in the neighbourhood of a'^''. In the case without exchange of the wave 
functions, 3?(6ij) as well as Q{bij) behave smoothly at a'^'^ while this is true only for 3?(6ij) 
in the case with exchange of the wave functions. In this case, Q'(%) jumps from a certain 
finite value y to —y at a". Since the Q{bij) of discrete states are zero, a jump in the '^{bij) 
can not appear in this case. The 3?(&jj), however, show a dependence on a which is very 
similar to that of resonance states with exchange of the wave functions (7^ < jf). 

In order to trace the influence of the branch point in the complex plane onto the mixing 
of discrete states, the differences 6 = — l^i^yiP and the values \bij\'^ are shown in Figs. 
§ and 1^ for different values 7^ from 7^ > 7" to 7^ = 0. Most interesting is the change 
of the value 6 from 1 to at 7?'". The relation = at 7 < 7" is the result 

from interference processes. It holds also at 7^ = 0, i.e. for discrete states. In this case, 
l&iiP = l&ijViP = 0.5 at a". 
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The values A and B characterizing the bi-orthogonahty of the two wave functions are 
shown in Fig. ^ for the same values of 7, as in Figs. ^ and ^. The A and B are similar for 
7j — 'jf^ = ±A as long as A is small. They approach A ^ 1 and B ^ for 7^ — 0. 

In Fig. 0, the energies Ei and mixing coefficients are shown for illustration for four 
discrete states with three neighboured avoided crossings as a function of a. In analogy to 
(0), the matrix is 

/ei(a) \ 
(4) ^ 62(0) 
63(0) 
V 64(0) / 

The mixing in the eigenfunctions of 7i which is caused by the avoided crossings remains, at 
high level density, at all values of the parameter a. It is the result of complicated interference 
processes. This can be seen best by comparing the two pictures with four interacting states 
(top and middle in Fig. |^) with those of only two interacting states (bottom of Fig. ^ and 
bottom right of Fig. Fig. |^ bottom shows the large region of the a values around a™ 
for which the two wave functions remain mixed: 1 and l&yyjp — > for a ^ with 

— 034! 3> \a^l — 034! ; i = 1,2. The avoided crossings between neighboured states do, 
therefore, not occur between states with pure wave functions and it is impossible to identify 
the unequivocally (Figs, fright top and middle). These avoided crossings are caused 
by branch points which overlap in the complex plane while the avoided crossings considered 
in Figs. to^ and Fig. ^bottom correspond to isolated hi a.nch points in the complex plane. 



/ UJ12 UJ13 LOu \ 

UJ2I 0^23 ^^24 

^31 ^32 UJ^4 

Vt^41 u;42 U;43 / 



(46) 



V. DISCUSSION OF THE RESULTS 

Most calculations represented in the present paper are performed for two states which 
cross or avoid crossing under the influence of an interaction u which is real. A general 
feature appearing in all the results is the repulsion of the levels in energy (except in the 
very neighbourhood of a™ when 7j > 7^'^'', i.e. F™(a,Ci;) < 0). This result follows analytically 
from the eigenvalue equation (pSj). It holds quite generally for real u as shown by means of 
the spectra of microwave cavities [§] and laser-induced continuum structures in atoms ||16|| . 
The level repulsion in energy is accompanied by an approaching of the lifetimes (widths) of 
the states. 

The sign of -F"(a, cu), Eq. (^), is decisive whether or not the states will be exchanged 
at the critical value a'^^ of the tuning parameter. When u is real and so small that F'^'^ < 
and the difference of the widths Fj — Fj^j is different from zero at a*^"" then the states will 
not be exchanged and the energy trajectories cross freely. If, however, > and Fj = Tj^i 
at a'^'^, the states will be exchanged and the energy trajectories avoid crossing. 

The exchange of the wave functions can be traced back to the branch point in the 
complex plane where the S matrix has a double pole and $j — ± i $j according to Eq. 
(^) (and calculations for a realistic case Ijl^)- Here, the real as well as the imaginary parts 
of the components of the wave functions increase up to an infinite value and (<l'*|$j) is the 
difference between two infintely large values. Thus, the orthogonality relation = 

and the normalization condition = 1 can not be distinguished. This makes possible 

the exchange of the two wave functions. 
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The exchange of the wave functions continues analytically into the function space of 
discrete states as illustrated in Figs. ^ § and ||. When the resonance states avoid crossing 
at a", the components of the wave functions do not increase up to infinity. Their increase is 
reduced due to interferences (Fig. |^). The differences S = |6mP— \hij^i\^ jump from 1 to at 
a" (and from about to almost 1 for values of a distant from a'^^) when {'^i — '^f^) changes its 
sign (Fig. ^). This jump is related to the exchange of the wave functions. The value 5 = 
at a'^'' remains unaltered when (71/2 — 72/2)^ < 4 cj^, i.e. also for discrete states. Therefore, 
the results shown in Fig. ^ top and middle correspond to situations being fundamentally 
and topologically different from one another. 

The bi-orthogonality ([22| ) of the wave functions is characteristic of the avoided crossing 
of resonance states. It increases limitless at the double pole of the S matrix and vanishes in 
the case of an avoided crossing of discrete states (Fig. ^). It does not enter any physically 
relevant values since it does not enter the S matrix, Eq. (0). The wave functions of the 
resonance states appear in the S matrix in accordance with the orthogonality relations (^Tj) 
which are fulfilled in the whole function space without any exceptions. 

Another result of the present study is the influence of the branch points in the complex 
plane onto the purity of the wave functions At a", the wave functions are not only 
exchanged but become mixed. The mixing occurs not only at the critical point a'^^ but in a 
certain region around a^"^ when the crossing is avoided. This fact is important at high level 
density where, as a rule, an avoided crossing with another level appears before $° is 

reached. As a result, all the wave functions of closely- lying states contain components of all 
basic states. That means, they are strongly mixed at high level density (for illustration see 
Fig. 0). 

The strong mixing of the wave functions of a quantum system at high level density means 
that the information on the individual properties of the discrete states described by the $° 
is lost. While the exchange of the wave functions itself is of no interest for a statistical 
consideration of the states, the accompanying mixing of the wave functions <I>i is decisive for 
the statistics. At high level density, the number of branch points is relatively large (although 
of measure zero). Therefore, the (discrete as well as resonance) states of quantum systems 
at high level density do not contain any information on the basic states with wave functions 
It follows further that the statistical properties of quantum systems at high level density 
are different from those at low level density. States at the border of the spectrum are almost 
not influenced by branch points in the complex plane since there are almost no states which 
could cross or avoid crossing with others states. The properties of these states are expected 
therefore to show more individual features than those at high level density. In other words, 
the information on the individual properties of the states with wave functions $° at the 
border of the spectrum is kept to a great deal in contrast to that on the states in the center 
of the spectrum. 

Thus, there is an influence of the continuum onto the properties of a (closed) quantum 
system with discrete states due to the analyticity of the wave functions. The branch points 
in the complex plane are hidden crossings, indeed. They play an important role not only in 
atoms, as supposed in [l^Jl^ , but determine the properties of all (closed and open) quantum 
systems at high level density. Their relation to Berry phases has been studied experimen- 
tally as well as theoretically [jll]-^ . 
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Summarizing the results obtained for quantum systems at high level density with avoided 
level crossings and double poles of the S matrix, it can be stated the following: 

- the poles of the S matrix correspond to the eigenvalues of a non-hermitian effective 
Hamilton operator also in the case that the resonance states overlap, 

- the eigenfunctions of a non-Hermitian Hamilton operator are bi-orthogonal in the 
whole function space without any exceptions, 

- avoided level crossings in the complex plane as well as in the function space of discrete 
states can be traced back to the existence of branch points in the complex plane, 

- under certain conditions, a branch point in the complex plane appears as a double 
pole of the S matrix, 

- branch points in the complex plane cause an exchange of the wave functions and create 
a mixing of the states of a quantum system at high level density even if the system is 
closed and the states are discrete. 

All the results obtained show the strong influence of the branch points in the complex 
plane on the dynamics of many-level quantum systems. They cause an avoided overlapping 
of resonance states which is accompanied by an exchange of the wave functions. A special 
case is the avoided level crossing of discrete states known for a long time. The avoided level 
crossings cause a mixing of the eigenfunctions of Ti which is the larger the higher the level 
density is. The states at the border of the spectrum of a many-particle system are therefore 
less influenced by avoided level crossings than those in the centre. 

Acknowledgment: Valuable discussions with E. Persson, J.M. Rost, P. Seba, M. Sieber, 
E.A. Solov'ev and H.J. Stockmann are gratefully acknowledged. 
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FIG. 1. The energies Ei (top left) and widths Ti/2 (top right) of the two eigenstates of the 



matrix (34) as a function of the parameter a. The thin hnes give the energies Ei and widths Ti/2 
of the states at a; = 0. The lower part of the figure shows the coefficients ba (bottom left) and bij^^i 
(bottom right) defined by Eq. (|37|). The x and o denote the ^{hj) while the 9(6ij) are denoted 
by + and o. ei = 1 — a/2; 62 = a; 7i/2 = 1.0; 72/2 = 1.1 and w = 0.05. 
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a a 

FIG. 2. The energies Ei (left) and widths Ti/2 (right) as a function of the tuning parameter a. 
ei = l-a/2; 62 = a and a; = 0.05. The 71/2 are 1.10 (top), 0.90 (middle), (bottom); 72 = 1.1 -71. 
The full lines show the and rj/2 for a; = 0. 
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a a 

FIG. 3. The mixing coefficients bu (left) and bij^i (right) defined by Eq. (|37| ) as a function 
of the tuning parameter a. o and x denote the real parts and o and + the imaginary parts, 
ei = 1 — a/2; 62 = a and u> = 0.05. The 71/2 are the same as in Fig. |2|: 1.10 (top), 0.90 (middle), 
(bottom); 72 = 1.1 -71. Note the different scales in the three cases. 
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a a 

FIG. 4. The differences 5 = \bii\'^ — \bij^i\'^ as a function of the tuning parameter 
ei = 1 - a/2; eg = a and w = 0.05. The 71/2 are 1.010 (top left), 0.990 (bottom left), 0. 
(top right), (bottom right); 72 = 1.1 • 71. Note the different scales in the different figures. 
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FIG. 5. The (full lines) and (dash-dotted lines) as a function of the tuning 

parameter a. ei = 1 — a/2; 62 = a and io = 0.05. The 71/2 are the same as in Fig. ^: 1.010 (top 
left), 0.990 (bottom left), 0.90 (top right), (bottom right); 72 = 1.1 • 71. Note the different scales 
in the different figures. 
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FIG. 6. The A (full lines) and B (dash-dotted lines) defined in Eq. ( p2[ ) as a function of the 
tuning parameter a. ei = 1 — a/2; 62 = a and to = 0.05. The 7i/2 are the same as in Fig. 1.010 
(top left), 0.990 (bottom left), 0.90 (top right), (bottom right); 72 = 1.1 • 71. Note the different 
scales in the different figures. 
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FIG. 7. The energies Ei (left) and mixing coefficients (right) of four discrete states (7^ = 
for i = 1,...,4) obtained from n^'^\ Eq. (|6|), as a function of the tuning parameter a. Top and 
middle: ei = 1 — a/3; 62 = 1 — 5a/12; 63 = 1 — a/2; = a; uj = 0.05 (top) and 0.1 (middle) 
for all non-diagonal matrix elements. Bottom: the same as above but ei = 1; 62 = 1-2; u; = for 
the coupling between the states i = 1,2 and j i^i, lo = 0.1 for the coupling between i = 3, 4 and 
j = 4,3. In this case, > \bij^i\^ (bottom right) as in Fig. |5| (bottom right). The dash-dotted 
lines (left) show Ei for a; = 0. The states i and j are exchanged at some values a. 
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